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CHAPTER  I 
INTRODUCTION 

Advanced  composite  materials  and  in  particular  Graphite/ 
Epoxy  are  finding  increasing  use  in  aerospace  vehicles  includ- 
ing rotary  wing  aircraft,  high  performance  fighters  and  the 
space  shuttle.'*'  Among  their  advantages  they  give  to  the 
designer  the  ability  to  taylor  the  elastic  and  aeroelastic 
characteristics  of  a structure  without  changing  its  physical 

dimensions.  Such  tayloring  is  now  being  considered  in  the 

• 2 
design  of  high  aspect  ratio  wings,  and  of  gas  turbine  com- 
pressor and  fan  blades. ^ In  some  cases  this  tayloring  takes 
the  form  of  shifting  one  or  several  of  the  natural  frequencies 
of  a structure.  In  other  cases  the  change  is  desired  in  a 
given  mode  shape  or  deflection. 

This  study  will  investigate  the  natural  mode  shapes  and 
frequencies  of  laminated  cantilever  plates  and  cylindrical 
shell  sections,  which  can  be  thought  of  as  being  representative 
of  fan  blades  and  low  aspect  ratio  wings  (Fig.  1) . 

The  purposes  of  this  study  are  to  develop  a simple  method  for 
understanding  the  vibratory  behavior  of  laminated  plates  which 
can  be  used  to  make  preliminary  design  estimates  and  which  can 
serve  as  a basis  for  nondimensionalizing  test  data.  The  second 
goal  is  to  model  the  laminated  plates  and  shell  more  exactly 
using  the  assumed  stress  hybrid  finite  shell  element  RS40,  and 
to  investigate  the  effects  of  transverse  shear  stiffness  on  the 


x 


natural  frequencies  using  this  element.  The  final  goal  is  to 
build  a set  of  laminated  plates  and  shells  and  test  them  for 
natural  frequency  and  mode  shape  to  check  the  accuracy  of  the 
finite  element  analysis. 

In  Chapter  II  the  free  vibration  of  a specially  orthotro- 
pic cantilevered  plate  is  analyzed  using  a single  assumed  mode 
partial  Ritz  analysis.  A form  of  the  solution  is  derived 
which  expresses  the  bending  and  torsional  frequencies  of  ortho- 
tropic plates  in  terms  of  plate  stiffnesses  and  known  values 
for  the  frequencies  of  isotropic  plates.  The  form  of  the 
solution  is  then  extended  to  chordwise  frequencies,  and  is 
( suggested  as  a basis  for  the  nondimensionalization  of  finite 

element  and  experimental  data. 

Chapter  III  discusses  the  finite  element  analysis  perform- 
ed using  two  seperate  elements;  an  assumed  stress  hybrid  shell 
element  RS40  and  a multilayered  plate  element  MLP3K. 

3 

The  software  used  to  implement  these  elements,  as  well  as  their 
convergence  properties  and  frequency  results  are  presented. 

The  frequencies  are  nondimensionalized  using  the  method  of 
Chapter  II,  and  correlated  to  the  isotropic  frequencies. 

The  selection  and  fabrication  of  the  test  articles  are 
presented  in  Chapter  IV,  along  with  the  experimental  techniques 
and  procedures  used.  Finally  in  Chapter  V the  experimental 
and  finite  element  results  are  compared  and  discussed. 
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CHAPTER  II 

Partial  Ritz  Analysis  for  Frequency  Determination 

2 . 1 Motivation 

Before  beginning  either  a detailed  finite  element  analysis 
or  an  actual  design  of  a laminated  cantilevered  plate  it  would 
be  desirable  to  have  some  feel  for  the  effect  of  ply  stacking 
sequence  and  orientation  on  the  mode  shapes  and  frequencies. 
Ideally  this  could  be  done  by  an  exact  solution  of  the  govern- 
ing partial  differential  equation,  but  since  the  exact  solution 
for  the  natural  modes  of  a specially  orthotropic  plate  have 
been  found  only  for  plates  with  2 or  4 sides  simply  supported, 
we  must  use  an  approximate  method.  In  this  case  a Partial 
Ritz  (also  called  Kantorovich  method) 4 technique  will  be  used 
which  assumes  a mode  shape  only  in  the  chordwise  direction, 
thus  reducing  the  governing  partial  differential  equation  to 
an  ordinary  differential  equation.  By  using  only  one  assumed 
mode,  the  approximate  form  and  functional  dependence  of  the 
frequencies  will  be  determined. 

2.2  Partial  Ritz  Analysis 

We  will  consider  the  transverse  vibration  of  a uniform 
specially  orthotropic  cantilever  plate  representative  of  a 
symmetric  laminated  plate  with  no  bending- twisting  coupling 
(Dig  = = 0) . Under  these  conditions  the  strain  energy  and 

kinetic  energy  are  given  as : ^ 
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V = 


\ ff  { 

* A 


_ .32w.2  , 32w 

D11  2 + 2D12  2 
J"L  3x  3x 


lf«  + D (if")  2 
. 2 22  \ 2 ' 
3y  3y 


+ 4D66(fw)2  1 “ 


(2.1) 


T = | //  ph(|^)2  dA  (2.2) 

A 

Applying  Hamilton's  principle  and  letting  the  harmonic  dis- 
placements of  free  vibration  take  the  form 


w(t)  = w sin  wt 


(2.3) 


we  find  the  following  equilibrium  condition  for  free  vibration: 


_ 1 f f r n , 32wN  2 ,3  w 3 w>  , n .3 w. 2 

2 * Dll(„  2)  + 2D12(7TT  77?  D22 


3x 


3x  3y 


3y 


+ 4D66 (gx^y) 2 “ W2phw2  } dA  = stationary  value  (2.4) 


Let  us  assume  a single  mode  in  the  chordwise  direction, 
so  that  the  displacement  everywhere  can  be  written: 


w(x,y)  = w(x)  + y 0(x) 


(2.5) 


where  Eq.  2.5  approximately  represents  the  bend  and  twist  of  a 
prismatic  bar.  Substituting  the  assumed  mode  shape  into  the 
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equilibrium  expression  (Eq.  214)  yields: 


Jt  - 

n = \ / f2  ( Di;l(w"  + yQ" ) 2 + 4D66(0')2 


2—  — 2 

- to  ph(w  + yO)  } dydx  = stationary  value 


(2.6) 


Note  that  the  stiffnesses  D22  and  D^2  have  been  eliminated 
because  no  curvature  is  allowed  in  the  chordwise  (y)  direction. 
Performing  the  integration  in  the  y direction,  it  is  convenient 
to  define  the  following  integrals,  the  limits  of  integration  of 
which  go  from  -c/2  to  c/2: 


a1(x)  = / D-^dy 

a2  (x)  = / Di;l  y dy 

2 

a3(x)  = / Di:l  y dy 

a4(x)  = f °66  dy 


oi2m^ (x) 
to2m2  (x) 
<j2m3  (x) 


/ tj2ph  dy 
2— 

f (u  ph  y dy 
/ w2ph  y2  dy 


(2.7) 


For  a plate  symmetric  about  the  mid  chord  point  (y  = 0) 
which  is  our  case,  a2  = m2  = 0 (Fig.  1) . In  this  case  the 
expression  for  the  potential  energy  reduces  to: 


a;L(w")  2 + a 3 ( 0 " ) 2 + 4a4(0')2 


2 -2 
ui  m^w 


2 n2 

- a)  m30 


dx  = stationary  value 


(2.8) 


Applying  the  variation  principle  to  Eq.  2.8  such  that 


6 it  = 0 


(2.9) 


Preforming  the  required  integration  by  parts,  and  gathering 
like  terms,  we  end  up  with: 


x.  _ 2-  - 

6u  = / { [a.w""  - m.u  w]  6w 

o 1 


+ [ a 3 0 " " - 4a40"  - u m30]  <$0  } dx 


Z Z 

+ [a.w"  ] 6w'  | - [a.w'"]  6w  | 

1 n x n 


Z Z 

+ [a~0"]  60'  | - [ a o 0 " 1 - 4a40']  60  | = 0 

■*  0 0 


(2.1Q) 


which  yield  the  governing  uncoupled  ordinary  differential 
equations  and  boundary  conditions  for  a plate  unsupported  along 
y = ±^.  In  the  case  of  a uniform  plate  clamped  at  x = 0 and 
free  at  x = Z,  the  appropriate  equations  and  boundary 
conditions  are: 


— 2—  — 
DiiCw""  - w phew  = 0 


J 2—  ~ J 

D. , 0""  - 4D,  ,-C0"  - u ph  © 


11  12 


(2.11) 


(2.12) 


( 
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2 .4 

d = UP.L  (2.20) 

11 

Eqs.  2.19  and  2.17  show  how  the  effect  of  the  warping 
stiffness  decreases  as  Dgg/D^j_  and  the  plate  aspect  ratio 
increase.  To  solve  Eq.  2.17  let 

0 = 0Qepx  (2.21) 

Then  the  characteristic  equation  becomes: 


( 


4 2 

p - b p - d = 0 

2 

a quadratic  in  p . Then  we  find 


2 _ b±/b2  + 4d 
P = 5 


(2.22) 


(2.23) 


Since  both  b and  d must  be  positive  numbers,  the  term 

2 

b + 4d  is  always  positive  and  greater  than  b,  and  the  values 
for  p can  be  given  as: 


( 


2 A>2  + 4d  + b c2 

p = ^ = f 

(2.24) 

_2  ,/b2  + 4d  - b,  2 

P = - ( 2 5 = ~ 9 

Then  the  solution  for  0{x)  can  be  written  in  either  of 


the  forms : 
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0(x)  = A1efi  + A2e-fx  + A3eigx  + A4e~igx  (2.25) 


or 


0(x)  = cos  gx  + B2  sin  gx 


+ B3  cosh  fx  + B^  sinh  fx 


(2.26) 


To  evaluate  the  frequencies  for  this  assumed  torsional 
mode,  one  substitutes  Eq.  2.26  into  the  cantilever  boundary 
conditions  on  0 given  in  Eq.  2.18,  and  the  result  is  a 4 by  4 
set  of  linear  homogeneous  equations  of  the  form: 


( 


[H]  {B}  = {0} 


(2.27) 


Or  written  explicitly: 


1 

0 

2 

-g  cosg 


0 

g 

2 . 

-g  sing 


1 

0 

f1 2coshf 


0 

f 

f 2sinhf 


(g  +bg)sing  - (g  +bg) cosg  (f  -bf)sinhf  (f  -bf)coshf 


*1 

" “ 

B1 

0 

Z. 

0 

B3 

0 

B4 

0 

(2 

28) 

For  a given  plate  the  value  for  b is  known  from  Eq.  2.19,  and 
f and  g are  related  by  Eq.  2.24  as: 


< 


~2  2 
f = g + b 


L. , 


(2.29) 
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Thus  given  a value  of  b,  g is  the  only  unknown  in  the  H matrix 
of  Eq.  2.28.  The  eigenvalues  of  Eq.  2.11  will  then  be  those 
values  of  g which  force  the  determinant  of  H to  zero.  The  re- 
sulting expressions  are  quite  lengthy  and  would  have  to  be 
solved  numerically  for  any  practical  problem.  Since  the 
assumed  mode  is  only  a rough  approximation  of  the  true  behavior 
of  the  plate,  the  accuracy  of  such  a result  would  probably  be 
less  than  that  predicted  by  a finite  element  or  finite  differ- 
ence model. 

However  it  is  very  instructive  to  look  at  the  form  that 
solutions  to  Eq.  2.28  will  take.  If  we  look  at  the  expression 
resulting  from  the  imposition  of  the  homogeneous  boundary 
conditions  at  x = 0 and  x = l in  this  problem,  we  get  express- 
ions of  the  form 


g = 


k 


n 


and  f 


+ b 


- £ 


n 


+ b 


(2.30) 


where  kn  is  the  eigenvalue.  Substituting  from  Eq.  2.24,  2.20 
and  2.19: 


,/43D66Ar2,2  . 4cAi4  ,48D66Ar2 

V D11  D11  Dil 


(2.31) 


which,  when  solved  for  to  gives  a relation  of  the  form 


(0 


2 


+ 48D66> 


(2.32) 
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Putting  the  expression  in  a form  similar  to  Eq.  2.16 


(2.33) 


(2.34) 


Thus  this  expression  indicates  that  the  torsional  fre- 
quencies of  a cantilever  plate  depend  linearly  on  the  aspect 
ratio,  and  upon  the  square  root  of  a "weighted  average"  of  the 
bending  (i.e.  warping)  stiffness  and  the  St.  Venant  torsional 
stiffness.  Further  the  influence  of  the  bending  stiffness  on 
the  weighted  average,  D,p,  increases  as  the  square  of  the  model 
eigenvalue,  and  decreases  as  the  square  of  the  plate  aspect 
ratio.  The  stiffness  does  not  appear  in  the  expression  for 

Dt  due  to  the  choice  of  a chordwise  mode  shape  which  excludes 
chordwise  bending.  If  a more  exact  solution  of  the  problem 
were  found,  one  might  expect  the  constants  in  Eq.  2.34  to 
change  and  perhaps  the  effects  of  to  be  included,  but  the 

general  form  of  Eq.  2.34  and  2.33  would  remain  unchanged. 

The  roughly  linear  dependence  of  wn  on  aspect  ratio  is 
well  born  out  by  more  complete  Ritz  analysis  and  by  experimen- 
tation for  isotropic  cantilever  plates.^  We  will  return  to 
the  exact  dependence  of  for  orthotropic  plates  after  the 
presentation  of  finite  element  analysis  results  in  Chapter  III. 
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2.4  Generalization  of  the  Frequency 
Nondimentionalization  Scheme 

Equations  2.16  and  2.33  strongly  suggest  that  the  natural 
frequencies  of  a cantilever  orthotropic  plate  can  be  correlated 
to  its  bending  stiffnesses  D.^ , D^2  / etc,  by  use  of  the  well 
known  expression  for  frequencies  of  free  vibration  of  a plate 


(2.35) 


( 


where  now  D is  interpreted  as  the  "appropriate"  bending  stiff- 
ness  for  the  mode  in  question.  We  have  seen  that  for  bending 
D = D i , the  plate  stiffness  in  the  bending  direction. 

For  torsion  however,  D = D , a weighted  average  of  the  bend- 

3.  i 

ing  and  torsion  stiffnesses. 

The  extension  of  this  scheme  to  the  simple  chordwise  modes 
seems  quite  obvious.  Here  we  might  expect: 


(2.36) 


The  "appropriate"  stiffness  D is  new  taken  as  the  bend- 

cl 

ing  stiffness  in  the  chordwise  direction,  and  kn  is  dependent 
on  the  aspect  ratio  of  the  plate. 

The  validity  of  these  assumptions  will  be  born  out  by 
their  ability  to  predict  the  actual  natural  frequencies  of  a 
given  laminated  plate,  or  those  predicted  by  finite  element 


€ 


methods.  We  have  however  at  least  the  basis  for  a scheme  for 
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the  rational  nondimentionalization  of  the  frequencies  for 
these  plates,  and  a starting  point  for  the  reduction  of  the 
data  from  the  finite  element  analysis  and  experiments. 


i 
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CHAPTER  III 

Finite  Element  Analysis 

3.1  Choice  of  Elements  and  Specimens  to  Model 

In  the  selection  of  a finite  element  to  use  in  modeling 
any  physical  system,  one  of  the  choices  which  must  be  made  is 
in  the  relative  sophistication  of  the  element  to  be  used. 

In  the  case  of  modeling  a variable  thickness  doubly  curved 
laminated  wing  or  blade,  it  may  be  advantageous  to  use  an 
element  which  includes  transverse  shear  effects  and  allows  for 
surface  curvature.  Elements  of  this  type  include  the  four 
node  thick-shell  quadrilateral  isoparametric  element  QUAD  4 

7 

developed  by  MacNeal,  and  the  six  node  thin-shell  triangular 
. 3 

isoparametric  element  developed  by  Minich.  In  this  study  the 
eight  node  quadrilateral  element  RS40  was  used  to  model 
graphite/epoxy  plates  and  shells.  To  serve  as  a basis  for 
comparison  the  same  structures  were  also  modeled  with  a plate 
element  KLP3K.  These  elements  will  be  discussed  in  detail  in 
the  next  section. 

The  choice  of  structures  to  model  was  largely  determined 
by  what  could  be  fabricated  and  tested,  and  will  be  discussed 
in  section  4.1.  For  reference  at  this  point  the  samples 
modeled  and  tested  included  3"  x 3"  flat  plates  of  G/E  and 
aluminum,  6"  x 3"  flat  plates  of  G/E,  G/E/Al,  and  6"  x 3" 
cylindrical  section  shells  of  G/E.  All  samples  were 
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cantilevered  (Fig.  1) . The  material  properties  used  in  the 
analysis  are  listed  in  Appendix  A. 

3.2  Description  of  Elements  MLP3K  and  RS4Q 

The  element  MLP3K  is  a 4 node  20  degrees  of  freedom 
constant  thickness  multilayered  plate  element.  It  is  available 
in  either  a triangular  or  quadrilateral  version,  and  is  based 
on  a moderately  thick  plate  formulation.  It  has  been  used  with 

g 

good  results  to  model  both  thin  and  moderately  thick  plates. 

In  this  study  the  quadrilateral  version  was  used  as  a rectang- 
ular element  to  model  thin  (t/c  ~ 1/70)  plates.  Element  MLP3K 
uses  linear  interpolation  of  displacements  on  its  boundaries, 
and  is  not  a compatible  element. 

The  element  stiffness  matrix  is  based  on  a hybrid  stress 
model  which  used  the  modified  complimentary  energy  (irmc^) 
variational  principle.  Since  the  element  models  a multilayer 
plate  a stress  field  is  assumed  within  each  layer  in  terms  of 
a standardized  set  of  stress  distribution  parameters. 

Interlayer  stress  compatibility  is  satisfied,  as  are  the  stress 
boundary  conditions  on  top  and  bottom  surfaces.  The  element 
behavior  is  defined  in  terms  of  displacements  u,  v,  w and 
rotations  0x  and  0y  at  each  node.  With  4 nodes  and  5 degrees 
of  freedom  (d.o.f.)  per  node,  there  are  20  element  level  d.o.f.. 

Transverse  shear  effects  are  treated  in  an  average  sense. 
Lines  perpendicular  to  the  mid  surface  remain  straight  but  not 
necessarily  perpendicular  to  the  deformed  midsurface.  Thus  Gx 


and  Qy  correspond  to  but  are  not  identical  to  ~ and 
Layer  to  layer  cross  section  warping  is  not  allowed,  but  over- 
all warping  is  permitted  based  on  the  average  transverse  shear 
properties . 

The  inputs  to  the  element  consist  of  the  ply  orientation, 
thickness  and  density,  as  well  as  inplane  and  transverse 
elastic  properties  for  each  ply.  The  element  is  assumed  to  be 
uniform  and  of  constant  thickness  and  density. 

The  ply  properties  are  assembled  to  form  laminate  proper- 
ties, and  the  element  stiffness  matrix  is  given  as  an  output. 
There  are  no  restrictions  on  the  laminate  considered,  so  that 
unbalanced  (D16»  D2g  ^ and  unsymmetric  (B  = 0)  laminates  can 
be  analyzed  using  MLP3K.  A corresponding  routine,  MLP3M,  uses 
bilinear  interpolation  to  form  a hybrid  rational  mass  matrix 
which  includes  the  effect  of  rotary  inertia.  Reference  8 
contains  a more  detailed  description  of  MLP3K. 

The  element  RS40  is  a moderately  thick  linear  elastic 
shallow  shell  element.  It  is  quadrilateral  having  4 corner  and 
4 midside  nodes,  and  uses  quadratic  interpolation  of  displace- 
ments on  its  boundaries..  It  can  model  arbitrary  shallow  com- 
plex curvatures  and  variations  in  thickness  over  the  element. 
The  element  has  been  used  to  model  thin  and  moderately  thick 

Q 

shell  structures. 

The  element  is  based  on  the  modified  Hellinger-Reissner 
Principle  (tt  ) . It  is  designed  to  model  a homogeneous  mater- 
ial so  a preprocessor  must  be  used  when  RS40  is  employed  to 
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model  composite  structures.  The  inputs  to  the  preprocessor 
are  the  ply  thickness,  density,  orientation  and  elastic  proper- 
ties. Using  a suitable  lamination  theory,  the  preprocessor 
calculates  the  elastic  matrices  A,  B and  D of  the  laminate, 
the  transverse  shear  matrices,  and  the  translational  and  ro- 
tary inertia  at  each  node.  RS40  has  been  modified  for  use  with 
laminated  materials  such  that  it  accepts  these  values  rather 
than  just  the  material  E matrix.  It  uses  the  A matrix  when 
calculating  the  strain  energy  associated  with  extension,  D with 
bending,  B with  bending  stretching  coupling,  etc.  In  this  way 
RS40  can  correctly  calculate  the  stiffness  matrix  of  a shell 
composed  on  unbalanced  (D^g,  D2  6^^  or  unsY'minetr;i-c 
laminate.  No  reduced  stiffness  approximation  must  be  applied 
to  account  for  the  bending-stretching  coupling. 

As  in  MLP3K,  transverse  shear  effects  are  handled  in  an 
average  sense,  with  perpendiculars  which  remain  straight  but 
are  free  to  twist.  The  nodal  degrees  of  freedom  are  then  u,  v, 
w in  the  global  x,  y,  z coordinates,  and  Ox  and  Oy  about  local 
x and  y axes  which  are  mutually  orthogonal  and  are  tangent  to 
the  shell  midsurface.  With  8 nodes,  there  are  40  element 
degrees  of  freedom. 

The  input  to  the  element  is  now  for  each  node:  the  coor- 
dinate of  the  node,  a set  of  3 direction  vectors  tangent  and 
normal  to  the  shell  midsurface,  the  laminate  elastic  and  iner- 
tia properties  from  the  preprocessor,  and  the  shell  thickness. 
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The  properties  and  thicknesses  are  interpolated  biquadratically 
over  the  element.  The  outputs  consist  of  rhe  element  stiffness 
matrix  and  hybrid  rational  mass  matrix  which  includes  the 
effects  of  rotary  inertia.  Note  however  that  in  its  present 
formulation  RS40  cannot  model  the  case  of  a mass  distribution 
unsymmetric  about  the  midplane.  For  a more  complete  discussion 
of  RS40,  see  Reference  9. 

3.3  Finite  Element  Software 

The  software  used  for  the  analysis  of  cantilever  plates 
and  cylindrical  shell  sections  using  the  MLP3K  and  RS40  ele- 
ments was  based  on  the  Finite  Element  Analysis  Basic  Library, 
version  5 (FEABL5)  developed  by  the  Aeroelastic  and  Structures 

O 

Research  Lab  at  M.I.T.  FEABL  is  a series  of  modular  sub- 
routines which  are  linked  together  by  the  user's  main  program. 

The  flow  diagram  used  in  conjunction  with  RS40  is  shown  in 
Fig.  2.  After  initialization  of  variables  and  arrays,  subrou- 
tine INPUT  reads  the  parameters  which  specify  the  artical  to  be 
modeled:  ply  orientation,  thickness,  stacking,  dimentions  etc; 

and  of  the  solution  required:  grid,  convergence  criteria,  etc. 
Routines  SETUP  and  ORK  are  FEABL  supplied  housekeeping  routines. 
LAMIN  and  CLT  prepare  the  laminate  properties  for  use  by  RS40 
according  to  classical  laminate  theory  as  described  in  Section 
3.2.  RS40  then  calculates  element  level  k and  m matrices,  and 
with  clamped  boundary  conditions  applied  along  y = 0,  ASEMSL 
creates  the  global  K and  M properties. 


The  resulting  eigenvalue  problem  is  solved  by  subroutine 
SSPACE  using  the  subspace  interation  method  (SIM) . ^ The  SIM 
technique  reduces  the  full  space  problem  to  a representative 
subspace,  solves  the  eigenvalue  problem  in  the  subspace,  then 
recreates  the  full  space  eigenvectors. 

For  a problem  with  N unconstrained  degrees  of  freedom, 
the  procedure  works  as  follows.  The  user  specifies  P,  the 
number  of  eigenvectors  to  be  included  in  the  subspace. 

The  value  of  P is  chosen  in  the  range 

p + 3 < P < 2p  (3.1) 

where  p is  the  number  of  eigenvectors  for  which  convergence  is 
desired.  The  program  then  uses  an  initial  guess  for  the  set 
of  P eigenvectors  to  reduce  the  full  space  K and  M matrices  to 

their  subspace  equivalents: 

« 

. T . T . 

k = U(i)  K LJ(l)  m = U{1)  N U(l)  (3.2 

where  U ^ is  the  N x P matrix  of  eigenvectors.  The  resulting 
auxiliary  P x P eigenvalue  problem 

k u = m u c (3.3 

is  solved  in  routine  SSPACE  using  an  iterative  Jacobi  diagonal- 
ization  scheme.  At  the  end  of  the  step,  the  full  space 
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eigenvectors  are  recomputed  according  to 


u(i  + 1)  _ K“1 


M U 


(i) 


u 


(3.4) 


and  the  procedure  is  begun  again  with  Eq.'s  3.2  using  the  new- 
ly  computed  vectors.  Convergence  of  a mode  occurs  when  for  a 
given  mode 


(i) 


A3 


(i  + 1) 


l/U- 


(i) 


< e 


(3.5) 


where  the  allowable  tolerance  parameter  e is  specified  by  the 
user.  Here  the  subspace  eigenvalue  is  identical  to  the  full 
space  value. 

Once  the  solution  has  converged,  SSPACE  prints  the  re- 
quested eigen  values  and  vectors  and  returns  command  to  main. 

An  M.I.T.  maintained  subroutine  CONTUR  was  used  to  contour 
plot  the  mode  shapes  using  a line  printer.  The  program  used 
in  conjunction  with  MLP3K  was  nearly  identical  to  the  described 
except  that  it  could  analyze  only  flat  plates.  Subroutine  CLT 
was  not  used  since  MLP3K  does  lamination  theory  calculations 
internally. 

3.4  Convergence  and  Accuracy 

Since  a study  of  convergence  and  accuracy  implies  that  we 


( 


know  an  exact  solution  to  which  a result  should  converge,  all 
convergence  studies  were  done  modeling  a flat  rectangular 
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aluminum  plate.  For  this  problem  the  assumed  mode  shape  Ritz 
solution's  listed  in  Tables  4.49  and  4.52  of  Reference  6 were 
considered  the  reference  solutions. 

In  studying  the  grid  convergence  of  a cantilever  plate 
model  there  are  two  parameters,  of  interest,  the  number  of 
elements  in  the  spanwise  (x)  direction  and  chordwise  (y) 
directions.  Convergence  of  each  frequency  is  most  easily  seen 
if  one  fixes  the  chordwise  number  of  elements  and  varies  the 
spanwise  number.  Families  of  convergence  curves  are  then  pro- 
duced. Figure  3 shows  two  typical  sets  of  curves,  for  the  2nd 
bending  and  2nd  torsion  modes  of  a plate  of  aspect  ratio  two. 
For  each  mode  the  convergence  for  MLP3K  (with  5 elements  in 
the  y direction)  and  RS40  (with  3 elements  in  the  y direction) 
are  plotted  versus  the  total  number  of  unconstrained  degrees 
of  freedom.  Since  RS40  has  8 nodes  per  element,  and  MLP3K  has 
only  4,  the  total  number  of  degrees  of  freedom  in  the  y direc- 
tion are  roughly  equivalent  in  these  two  cases. 

Figure  3 demonstrates  that  RS40  quickly  converges  to  a 
value  between  1 and  2 percent  below  the  Ritz  solution.  On  the 
other  hand  MLP3K  converges  much  more  slowly,  probably  due  to 
the  linear,  rather  than  quadratic  interpolation  of  displace- 
ments used.  In  the  second  torsion  mode  MLP3K  converges  to  a 
frequency  about  4%  less  than  the  Ritz  value,  and  thus  repre- 
sents the  structure  as  being  too  soft.  The  plots  in  Fig.  3 are 
very  representative  of  the  other  modes,  and  for  plates  of 


AR  = 1.  In  almost  all  cases  RS40  converged  more  quickly  to  a 
more  accurate  value  than  did  MLP3K. 

The  final  choices  for  the  grids  used  represent  a tradeoff 
between  computing  accuracy  and  cost,  and  are  summarized  below: 


PLATE 


ELEMENT  ' 
MLP3K 
RS40 


3"  x 3" 

6”  x 3" 

7x5 

10  x 5 

3x3 

6x3 

Based  on  their  agreement  with  the  Ritz  analysis  for  an  isotro- 
pic plate  as  shown  in  Table  1,  the  predicted  values  of  fre- 
quency from  MLP3K  should  be  within  +3%  and  from  RS40  within 
+1%.  These  grid  choices  were  then  extended  to  all  laminate 
samples  modeled.  Isotropic  cylindrical  shell  sections  were 
also  modeled  with  RS40.  The  results  were  found  to  be  in  good 
agreement  with  simple  beam  theory  for  first  bending  modes,  and 
with  torsion  theory  for  torsional  modes. 

The  only  other  choice  left  to  the  user  is  the  number  of 
modes  for  which  convergence  is  desired,  and  the  number  of  modes 
to  be  included  in  the  subspace.  Since  convergence  was  desired 
for  5 modes,  runs  were  made  with  S,  9,  and  10  in  the  subspace 
with  no  difference  in  results.  In  all  remaining  calculations  8 
modes  were  kept  in  the  subspace. 

For  many  choices  of  grid  size  it  became  obvious  that  the 
subspace  iteration  produced  incorrect  values,  high  and  low,  for 
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the  first  bending  frequency.  This  problem  occurred  more  often 
when  the  ratio  of  the  number  of  chordwise  elements  to  that  of 
elements  was  larger  than  1/2. 

The  problem  was  also  sensitive  to  element  aspect  ratio. 

It  is  suspected  that  this  error  is  due  to  a bad  choice  for  the 
initial  guess  at  the  eigenvectors  made  in  subroutine  SSPACE. 
Frequency  values  which  are  thought  to  be  incorrect  are  entered 
in  the  Tables  in  parenthesis. 

From  this  point  on,  all  results  discussed  will  be  those 
from  RS40.  These  are  largely  substantiated  by  the  MLP3K 
results,  but  the  RS40  results  should  be  more  accurate. 


3.5  Nondimentionalized  Frequency  Results 

The  computed  frequencies  from  RS40  are  shown  in 
Tables  2-8,  and  the  computed  mode  shapes  as  shown  in  Figures 
5 to  74.  The  mode  shapes  are  shown  as  contour  plots  with  node 
lines  shown  by  dark  lines  and  contours  of  constant  deflection 
shown  by  alternating  bands  of  white  and  grey  (characters) . 

The  clamped  edge  is  shown  at  the  bottom.  The  actual  numerical 
frequencies  will  be  discussed  in  more  detail  when  they  are 
compared  with  measured  values  in  Chapter  5. 

The  scheme  of  "rational  nondimentionalization" , that  is 
nondimentionalization  of  the  frequency  by  the  appropriate 
stiffness  terms  was  proposed  in  Section  2.4.  This  scheme  was 


i 


applied  to  the  frequencies  calculated  by  RS40. 


Iii  Tables  2 and 
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3 the  frequencies  and  rationally  nondimensionalized  frequencies 

are  listed  for  a series  of  laminated  plates  of  AR  = 1 and  2. 

After  seme  trial  and  error,  the  scheme  which  gave  the  best 

results  for  k was  as  follows.  Let  to  be  the  predicted 
n n 

frequency,  then  kn  the  nondimensionalized  frequency  is  given  by: 

k Pd~ 

“n  = 7?  V or  k 


(3.6) 


where 

i = length  of  the  plate 
p = mass  per  unit  area 

Da  = "appropriate"  stiffness  for  the  mode 
analyzed. 

The  appropriate  stiffneses  were  found  as 

(3.7) 

(3.8) 

(3.9) 

The  form  of  Eq.  3.7  and  3.8  are  just  those  found  in 
Eq.  2.16  and  suggested  in  Eq.  2.36.  The  form  of  Eq.  3.9  is 
quite  similar  to  Eq.  2.34  found  from  the  Ritz  analysis  and 


Da  " D11 


for  spanwise  bending: 
for  chordwise  bending:  D = D_0 

3L  » 12. 

for  torsion: 


Da  * DT  “ (DllH  + 2d12  + 4D66> 
AR 


,n 


+ 2) 
AR 
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repeated  here: 


k * 

D n ^ + 48D_ 
11  AR2  66 


(3.10) 


Both  are  weighted  averages  of  the  torsion  and  warping 

stiffnesses.  The  term  48Dgg  has  been  replaced  by  the  form  of 

the  coefficients  from  the  plate  bending  equation,  and  the 
2 

factor  of  kn  has  simply  been  replaced  by  the  torsional  mode 
number  n (n  = 1 for  1st  torsion,  2 for  second  torsion,  etc). 

The  term  in  the  denominator  of  Eq.  3.9  is  chosen  so  that  D 

cl 

goes  to  D for  an  isotropic  plate. 

The  frequencies  reduced  by  this  scheme  and  shown  in  Tables 
2 and  3 of  very  different  laminates  show  remarkable  consistancy 
and  agreement  with  the  isotropic  value.  In  every  case  listed 
in  Tables  2 and  3 if  one  estimated  the  frequencies  using 
Eq.  3.6  and  the  known  isotropic  value  for  k , the  result  would 
correctly  predict  the  ranking  of  the  frequencies  (1st  bending, 
1st  torsion,  etc). 

The  implication  of  this  is  that  using  the  simple  relations 
of  Eq.  3.6  through  3.9  and  classical  lamination  theory,  one  can 
estimate  the  relative  frequencies  and  approximate  mode  shape  of 
almost  any  thin  laminated  cantilever  plate  to  within  about  104. 
Since  the  lower  frequencies  are  usually  v/ell  spread  out,  this 
will  usually  predict  the  exact  ranking  of  the  natural  frequen- 
cies . 


CHAPTER  IV 


Experimental  Equipment  and  Procedures 

4.1  Choice  of  Specimens  to  be  Tested 

In  choosing  what  types  of  specimens  to  construct  and  test 
for  comparison  with  finite  element  predictions,  there  were  four 
main  parameters  considered:  the  effect  of  aspect  ratio,  of 
curvature,  of  the  anisotropy  of  the  laminate,  and  the  effect  of 
a sandwich  or  hybrid  laminate. 

To  study  the  effect  of  aspect  ratio  on  the  accuracy  of 
computed  frequency  and  mode  shapes,  it  was  decided  to  build 
Graphite/Epoxy  (G/E)  plates  of  aspect  ratio  1 and  2.  The  widest 
specimen  that  could  be  tested  was  3 inches,  so  the  sides  were 
set  at  3"  x 3"  and  6"  x 3".  To  verify  the  ability  of  RS40  to 
model  curved  shells  a set  of  cylindrical  shell  section,  6"  x 3" 
with  a radius  of  curvature  of  5"  were  tested.  In  many  practical 
problems,  Graphite/Epoxy  is  not  used  alone,  but  is  bonded  to  a 
substructure.  To  simulate  this  case  a set  of  6"  x 3"  Graphite/ 
Epoxy/Aluminum  symmetric  sandwich  plates  were  built  and  tested. 

In  each  of  the  above  cases,  the  effects  of  the  choice  of 
laminate  is  of  interest.  In  the  limit  one  can  have  a laminate 
very  stiff  in  bending,  stiff  in  torsion  or  quasi-isotropic. 

After  preliminary  design  and  testing  the  following  8 ply 
laminates  were  chosen  to  represent  the  three  cases  of  interest: 
bending  stiff  [02/+30]g,  torsion  stiff  [+45/+45]g,  and  quasi- 
isotropic [0/+45/90]g.  The  samples  tested  can  be  represented 
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by  the  matrix  below,  where  the  number  indicates  the  sample 
aspect  ratio: 


G/E  plate 

[02/±30]s 

1,2 

[+/ +45 ] s 

1,2 

[0/+45/90] 

1,2 

G/E  shell 

2 

2 

2 

G/E/Al  plate 

2 

2 

2 

( 


In  the  case  of  the  aluminum  sandwich,  the  8 plys  were 
placed  symmetrically  about  the  aluminum,  4 on  a side,  and  0^ 
replaced  the  02/+30. 

To  serve  as  a calibration  for  the  test  procedure  a set  of 
aluminum  plates  of  aspect  ratio  1 and  2 and  roughly  the  same 
thickness  were  built  and  tested. 

4.2  Fabrication  of  the  Specimens 

The  materials  used  in  fabricating  the  specimens  were  pri- 
marily Graphite/Epoxy  prepreg  and  rolled  aluminum  sheet. 

The  graphite  was  Hercules  type  AS/3501-6  taken  from  lot  676. 

The  aluminum  was  2024-T3  rolled  sheet  nominally  .040"  thick. 

The  lay  up  and  cure  process  depended  on  the  type  of  sample 
being  made.  For  the  flat  G/E  plates,  the  prepreg  was  layed 

•up  as  laminates  roughly  12"  x 14"  in  size.  The  laminate  was 
« 

cured  in  a press  using  the  cycle  listed  in  Appendix  B,  and 
surrounded  by  a symmetric  arrangement  of  peel  ply  and  bleeders. 
In  the  press  cycle,  the  manuf actureres  temperature  and  pressure 
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cycle  are  approximated,  but  no  vacuum  was  drawn  over  the 
laminate.  After  curing,  specimens  were  cut  from  the  laminate 
using  a diamond  coated  cutting  wheel. 

For  the  Graphite/Epoxy /Aluminum  sandwich  specimens,  the 
prepreg  was  cut  to  the  exact  size  of  the  aluminum  plate,  the 
plate  being  used  as  a template.  Six  of  these  6"  x 3"  samples 
were  then  cured  at  a time  in  the  press  with  the  graphite 
cocured  onto  the  aluminum.  The  samples  were  seperated  during 
curing  by  resin  dams  which  insured  the  resin  would  flow  verti- 
cally and  not  sideways,  and  that  the  fibers  would  not  flow 
arround  the  edges  of  the  plate.  The  same  cure  cycle  was  used 

as  with  the  flat  G/E  laminates. 

( . . 

The  cylindrical  shell  section  specimens  were  laved  up  in- 
dividually on  a cylindrical  tool.  The  lay  up  was  directly  on 
the  surface  of  the  tool.  The  laminate  was  surrounded  by  a 
resin  dam,  and  peel  ply  and  bleeders  were  stacked  on  top  of  the 
laminate  only.  The  curved  specimens  we re  cured  in  an  oven 
under  vacuum  using  the  second  cycle  in  Appendix  B.  This  cycle 
approximates  the  temperature  and  vacuum  aspects  of-  the  manu- 
facturers, but  the  pressure  on  the  laminate  is  only  about 
14  psi,  instead  of  the  100  psi  that  is  called  for. 

The  next  step  in  the  preparation  of  all  the  samples  was  to 
finish  the  edges  by  sanding.  The  edges  were  sanded  until  they 
were  parallel,  square  and  of  the  correct  dimensions.  The  final 
f step  before  the  sample  could  be  tested  was  to  bond  onto  the  base 
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a 1"  x 3"  loading  tab.  This  tab  was  to  insure  that  the  speci- 
men could  be  clamped  firmly  and  squarely  in  the  test  jig  without 
any  damage  being  done  to  the  surface  fibers.  The  aluminum  tabs 
were  machined  to  size  (and  to  the  proper  curvature  for  the 
curved  shell  specimens)  and  bonded  to  the  laminates  in  a 
secondary  bond  using  Cyanamid  FM123  film  adhesive  at  225 °F  and 
50  psi. 

Because  this  bonding  operation  introduced  a possible 
source  of  error  into  the  test,  half  of  the  aluminum  samples 
made  had  tabs  bonded  on  in  a manner  identical  to  those  on  the 
laminated  specimens.  The  other  half  were  machined  from  a solid 
plate  of  3/8"  aluminum  to  the  same  dimensions  as  the  bonded 
samples.  Since  no  difference  was  found  in  the  frequencies  and 
mode  shapes  of  the  bonded  and  integral  specimens,  the  bonded 
tab  is  assumed  to  introduce  no  noticible  error  in  the  testing. 

As  another  check  some  of  the  samples  were  made  symmetric  about 
the  clamp,  while  others  extended  only  in  one  direction  from  the 
clamp.  Again  no  significant  difference  was  found,  so  the 
asymmetry  of  the  specimens  about  the  clamp  is  assumed  to  have 
no  adverse  effect  or  the  results. 

4.3  Test  Apparatus  and  Proceedures 

The  test  apparatus  is  shown  in  Fig.  4.  The  specimen  was 
clamped  into  a rectangular  block  of  aluminum  1"  x 6"  x 6"  by 
three  3/8"  bolts  tightened  to  a torque  of  75  in#.  The  block 
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was  suspended  by  4 spring  steel  flexures.  The  flexures  allowed 
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the  block  to  move  perpendicular  to  the  clamp  but  were  very 
stiff  in  other  directions.  The  block  was  rigidly  attached  to 
a Ling  model  420  shaker  mounted  horizontally.  The  shaker  has 
a peak  force  output  of  100  lb.  force.  With  this  setup  the 
samples  could  be  clamped  very  rigidly  at  their  base  and  ex- 
cited along  their  z axis  in  a range  from  about  5 to  3000  Hz. 

The  moving  base  was  instrumented  with  2 Endevco  model 
7701-50  "Isoshear"  accelerometers.  One  Endevco  2222B  micro- 
ministure  accelerometer  was  mounted  on  the  specimen  within 
3/4"  of  the  base  and  centerline.  The  presence  of  this  sample 
mounted  accelerometer  was  found  to  have  less  than  a .2%  effect 
on  the  frequency,  and  often  no  effect  at  all. 

By  monitoring  the  base  and  sample  mounted  accelerometers 
on  a scope,  one  gets  a direct  indication  of  the  amplitude  and 
phase  response  of  the  sample.  Resonances  were  taken  as  the 
point  of  90°  phase  shift  between  the  input  and  response 
accelerometers.  In  the  case  of  some  less  strongly  excited 
modes  where  a full  90°  phase  shift  did  not  occur,  the  combina- 
tion of  maximym  phase  shift  and  amplitude  was  used  to  determine 
the  resonance.  The  frequency  data  was  read  from  a digital 
counter  and  found  to  be  reproducable  to  about  .2  Hz  for  bending 
modes,  .5  Hz  for  chordwise  modes  and  1 Hz  for  torsion  modes, 
which  were  only  lightly  excited. 

The  nodal  pattern  of  a given  mode  was  found  by  slowly 


scanning  the  surface  with  a pointed  pencil  eraser,  while 
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watching  the  phase  and  amplitude  response  on  the  scope.  Nodal 
points  could  be  easily  identified  as  those  points  which  when 
touched  had  no  effect  on  the  phase  or  amplitude  response. 

These  points  were  recorded  on  tracing  paper,  and  the  resulting 
modal  patterns  are  easily  identifiable.  In  this  way  the  fre- 
quency and  mode  shapes  for  the  first  5 modes  of  all  the  speci- 
mens were  observed  and  recorded. 
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CHAPTER  V 

Comparison  of  Experimental  and  Finite  Element  Results 


5.1  Frequency  Results  for  Aluminum  Plates 

A total  of  eight  aluminum  plates  were  tested,  four  3"  x 3" 
and  four  6"  x 3".  Since  the  plates  varied  in  thickness  from 
.04"  to  .05",  the  measured  frequencies  had  to  be  corrected  for 
this  effect.  Fom  the  expression  for  the  frequency  of  an  iso- 
tropic plate 


0) 


n 


kn  /~D  _ kn  / E t3 
a2  i p i2  1/12  (i  + v2)pt 


12 (1+v) 


2- 

P 


} 


(5.1) 


we  see  the  frequency  depends  linearly  on  the  thickness.  Using 
this  relation  the  measured  frequencies  were  normalized  to  a 
reference  thickness  of  .0416",  and  averaged  for  each  mode. 

The  resulting  values  are  shown  in  Table  4.  The  coefficient  of 
variation  of  nominally  identical  specimens  averaged  about  1% 
and  did  not  exceed  2.5%  for  any  mode,  showing  that  the  results 
were  fairly  consistent. 

The  frequency  results  predicted  by  RS40  finite  element 
analysis  for  an  aluminum  plate  .0416"  thick  are  also  listed  in 
Table  4.  In  comparing  the  measured  and  calculated  frequencies 
we  see  that  at  most  the  two  differ  by  2.0%  of  the  measured 
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frequency.  Furthermore  six  of  the  ten  modes  measured  for  the 
two  sizes  of  plates  are  within  0.9%.  This  represents  quite 
excellent  agreement  between  calculated  and  measured  values. 

With  this  good  agreement  for  isotropic  plates,  we  have 
some  confidence  that  the  frequencies  measured  for  the  laminated 
specimens  will  truely  represent  the  cantilever  frequencies,  and 
will  not  be  affected  by  any  imperfection  of  the  clamp,  or  by  the 
presence  of  the  bonded  loading  tab. 

5.2  Frequency  Results  for  the  Laminated  Plates  and  Shells 

Thirty-seven  laminated  plates  and  shells  were  tested,  an 
average  of  three  for  each  configuration  described  in  Section 
4.1.  In  no  case  were  more  than  5 nor  less  than  2 of  a given 
type  tested.  The  8 ply  Graphite/Epoxy  laminates  ranged  from 
.041"  to  .043"  in  thickness.  If  the  manufacturer  stated 
nominal  ply  thickness  of  .0052"  is  assumed,  the  8 ply  thickness 
should  be  .0416".  Since  this  was  the  value  used  in  the  finite 
element  analysis,  the  experimental  results  should  also  be 
normalized  to  this  thickness.  However  a correction  of  the  data 
to  account  for  fiber  volume  effects  implies  a knowledge  of  all 
the  ply  elastic  properties,  and  an  understanding  of  how  the  ply 
stiffness  contribute  the  effective  stiffness  of  a given  mode. 
Since  neither  of  these  are  well  known,  the  frequency  data  was 
corrected  for  thickness  using  the  simple  linear  relationship 
for  isotropic  plates  given  by  Eq.  5.1. 
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The  corrected  frequencies  were  averaged  and  are  listed  in 
Tables  5 to  8 for  the  various  laminates  and  configurations. 

The  results  were  remarkably  consistant,  with  a coefficient  of 
variation  which  averaged  1.0%  and  was  less  than  or  equal  to 
0.6%  for  27  of  the  60  modes  measured.  This  consistancy  of 
results  give  some  confidence  in  the  workmanship  and  quality  of 
specimens  tested. 

If  the  ranking  of  the  frequencies  given  in  Tables  5 to  8 
is  studied,  some  interesting  trends  can  be  observed.  Looking 
at  the  3"  x 3"  G/E  plates  in  Table  5,  we  note  that  the  "torsion 
stiff"  laminate  [+45/+45]_  most  closely  simulates  the  isotropic 

""  b 

case  in  spacing  and  ranking  of  frequencies,  while  the  laminate 
[0/+45/90]  which  is  "quasi  isotropic"  with  respect  to  in  plane 
properties  now  has  a different  ranking  than  the  aluminum,  with 
first  chordwise  lower  than  second  bending.  The  case  is  differ- 
ent in  the  6"  x 3"  plates  of  Table  6.  Here  the  [0/+45/90]s 
resembles  the  isotropic  case,  while  [+45/+45]g  has  become  flex- 
ible enough  in  bending  that  second  bending  is  now  lower  than 
first  torsion.  | 

As  shown  in  Table  7,  the  effect  of  the  aluminum  core  in 
the  sandwich  specimens  is  to  moderate  the  selective  stiffening 
effect  of  the  Graphite/Epcxy , and  to  decrease  the  anisotropy  of 
the  plate.  Dramatic  shifts  in  frequency  and  mode  shape  are  now 
harder  to  achieve.  The  [0/+45/90/A1] g laminate  now  has  indeed 
become  almost  isotropic  in  behavior.  The  introduction  of  only 


ilk 


( 


4.6%  curvature  (rise  height/chord)  has  a dramatic  effect  on 
the  frequencies  (Table  8) . The  slight  amount  of  curvature 
stiffens  the  bending  modes  greatly,  but  has  only  a small  effect 
on  the  torsion  and  chordwise  frequencies.  First  torsion  has 
now  replaced  first  bending  as  the  lowest  mode  for  the  curved 
samples  of  laminates  l0/+45/90]g  and  [02/+30]s,  and  second 
bending  has  disappeared  from  the  lowest  five  modes.  The  ability 
for  composites  to  taylor  the  modal  ranking  and  nodal  structure 
of  a cantilever  with  some  curvature  seems  even  greater  than  the 
effect  on  flat  plates. 

The  calculated  frequencies  shown  in  Tables  5 to  8 show 
fair  agreement  with  the  measured  values.  Some  values  are  with- 
in a few  percent  while  others  are  off  as  much  as  12%,  the 
average  error  being  about  6%.  In  every  case  the  ranking  of  the 
modes  was  correct,  and  except  for  1 value  in  60,  the  calculated 
values  are  too  high.  These  two  facts  suggest  the  basic  material 
elastic  constants  used  for  in  plane  and/cr  transverse  proper- 
ties may  be  too  high.  The  material  properties  used  and  listed 
in  Appendix  A are  the  design  values  used  by  the  Grumman 
Corporation  for  this  fiber/resin  system.  During  the  course  of 
this  study  other  investigators  at  M.I.T.  using  the  same  system 
and  cure  cycles  measured  material  constants  significantly 
different  than  those  listed  by  Grumman.  This  is  possible  due 
to  the  inability  to  reproduce  the  manufacturer's  recommended 
cure  cycle. 
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The  resulting  uncertainty  in  the  actual  material  pro- 
perties of  the  specimens  used  is  as  great  as  the  errors  in  fre- 
quency shown  in  Tables  5 to  8.  A more  realistic  set  of  material 
properties  was  used  to  recalculate  the  frequencies  in  a few 
cases,  as  shown  in  Table  9.  Here  the  second  set  of  properties 
in  Appendix  A was  used,  which  correspond  more  closely  to  those 
measured  at  M.I.T.  The  resulting  frequencies  agree  more  closely 
with  the  observed  values.  Thus  the  accuracy  of  the  RS40  finite 
element  calculations  is  difficult  to  evaluate  due  to  the  large 
uncertainty  in  material  properties. 

I 

During  an  attempt  to  achieve  a better  prediction  of  the 
observed  frequencies,  the  transverse  shear  stiffnesses  were 
among  those  varied.  It  was  observed  that  a 10%  reduction  in 
the  transverse  shear  stiffness  had  the  effect  of  reducing  the 
frequencies  much  less  than  1%.  It  seems  that  with  the  present 
treatment  the  softening  effect  of  transverse  shear  for  a thin 
cantilevered  plate  is  small. 

5.3  Mode  Shape  Results 

The  mode  shapes  as  experimentally  determined  are  shown  in 
Figs.  5 to  74.  The  clamped  edge  is  shown  at  the  bottom,  and 
the  solid  lines  indicate  the  node  lines,  (i.e.  the  points  of 
zero  displacement).  The  patterns  are  exactly  as  was  observed, 
and  have  not  been  corrected  in  any  way.  The  patterns  are  also 
labeled  with  the  number  in  assending  order  of  frequency  (1st  is 
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lowest,  5th  the  highest  measured)  and  their  common  name 
(1st  bending,  etc) . 

As  can  be  observed  the  agreement  in  the  calculated  and 
actual  nodal  pattern  is  excellent  in  almost  all  samples  and 
modes  observed. 
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CHAPTER  VI 
Conclusions 

6.1  Finite  Element  Analysis 

The  finite  element  RS40  is  an  effective  and  versatile 
anisotropic  shell  element  for  use  in  dynamic  modeling.  When 
used  with  a preprocessor  which  does  laminate  calculation,  RS40 
can  be  used  to  model  laminated  plates  and  shells  of  variable 
thickness,  mass  and  stiffness,  and  with  arbitrary  shallow  com- 
plex curvature.  The  element  has  good  convergence  properties, 
and  in  the  isotropic  case  agrees  with  analytic  and  experimental 
results  for  cantilever  plates. 

The  natural  f requencies^of  cantilever  composite  plates  and 
shells  predicted  by  element  RS40  were  in  reasonable  agreement 
with  the  experimentally  observed  values.  In  all  cases  the 
ranking  of  the  frequencies  was  correct.  In  view  of  the  uncer- 
tainty in  the  laminate  stiffness  properties,  it  is  difficult 
to  assess  the  accuracy  of  RS40.  The  predicted  mode  shapes 
are  in  excellent  agreement  with  the  observed  patterns. 

6.2  Nondimensional  Analysis 

The  approximate  values  and  rankings  of  natural  frequencies 
of  laminated  cantilevered  plates  can  be  determined  from  the 
known  values  of  for  an  isotropic  plate,  and  the  simple 
relation  ojn  =(k  /Jt^)  / Da/p  , where  is  the  appropriate  ortho- 
tropic plate  stiffness  for  the  mode  in  question.  Using  this 
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relation  and  classical  lamination  theory,  one  can  make  preli- 
minary design  judgements  and  evaluate  the  effect  of  ply  orien- 
tation and  stacking  sequence  on  natural  frequencies  without 
sophisticated  numerical  models. 

6.3  Use  of  Composites  for  Elastic  Tayloring 
Graphite/Epoxy  composites  can  be  used  effectively  to  alter 

the  natural  frequencies  of  a plate  or  shell  of  fixed  dimension, 
and  to  taylor  the  mode  shaper  to  produce  desirable  elastic  or 
aeroelastic  properties. 

6.4  Effect  of  Transverse  Shear 

Using  RS40,  the  softening  effect  of  transverse  shear  on 
the  frequencies  of  a thin  laminated  plate  were  found  to  be 
small. 


( 


49 


( 


APPENDIX  A 
Material  Properties 

Hercules  reported  data  for  AS/3501-6 
0°  Tensile  Modulus  20.0  x 10^  psi 

Design  properties  used  in  finite  element  calculations: 

(Based  on  Grumman  design  values.) 

El  18.5  x 106  psi 

6 


LT 


LT 


’13 


’23 


1.60  x 10 
.250 

.650  x 10* 
.650  x 10( 
.222  x 10* 


Modified  properties  used  in  calculating  frequencies  in 

Table  9,  and  based  on  M.I.T.  measurements: 

6 psi 
6 


el 

16.0 

X 

10 

et 

1.36 

X 

10 

VLT 

.275 

glt 

.900 

X 

10 

G13 

.650 

X 

10 

G2  3 

. 222 

X 

10 
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APPENDIX  B 
Cure  Cycles 

Cure  cycle  used  with  AS/3501-6  cured  in  hot  press: 


1. 

Preheat  press  to  275°F. 

2. 

Place  cold  workpiece  on  press. 

Apply  15  psi  pressure. 

3. 

Leave  at  275°F,  15  psi 

16 

min. 

(plates  heat  up  to  275°F) . 

4. 

Raise  temperautre  to  300 °F. 

5 

min. 

5. 

When  temperature  increases  to  300°F 

increase  pressure  to  100  psi. 

6. 

Leave  at  300°F,  100  psi. 

30 

min. 

7. 

Raise  temperaure  to  350°F. 

10 

min. 

8. 

Leave  at  350°F,  100  psi. 

30 

min . 

Cure  cycle  used  with  AS/3501-6  cured  in  oven: 


1. 

Place 

cold  workpiece  in  oven. 

2. 

Pull 

vacuum  to  14 

psi. 

3. 

Raise 

temperature 

to  275°F. 

25 

min . 

4. 

Leave 

at  275°F. 

30 

min. 

5. 

Raise 

temperature 

to  350 °F. 

15 

min. 

6. 

Leave 

at  350°F. 

60 

min . 

7. 

Turn 

off  power,  cool  down  in  oven. 

30 

min . 

8. 

Turn 

off  vacuum. 

remove  from  oven. 
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APPENDIX  C 

Formulation  of  the  Shell  Element  RS40 
The  shell  element  used  in  the  present  work  has  eight  nodes. 
Each  node  has  five  nodal  degrees  of  freedom.  The  element  is 
derived  from  a modified  Hellinger-Reissner  principle  which 
states  that"*" 
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where 


e:  inplane  strains  expressed  in  terms  of  displacement 
y:  transverse  shear  strains  expressed  in  terms  of 
displacement 

y:  independent  transverse  shear  strains 

C yi  elastic  constant  matrices 

W:  potential  energy  due  to  applied  load 

m:  mass  per  unit  volume 

u:  displacement  vector 

ii:  acceleration  vector 

V : volume  of  nth  element 
n 


The  finite  element  is  formulated  by  assuming  displacements 
in  terms  of  nodal  displacements  and  by  assuming  independent 
transverse  shear  strains  in  terms  of  unkown 


: 
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parameters.  The  assumed  strain  parameters  are  to  be  elimimited’ 
at  the  element  level.  For  the  description  of  geometry  and  dis- 
placement we  define  two  cartesian  coordinate  systems,  the  glo- 
bal coordinates  and  local  coordinates.  See  Fig.  c-1:  One  axis 
of  the  local  cartesian  coordinates  is  normal  to  the  shell  mid- 
surface  while  the  other  two  axes  are  tangent  to  the  midsurface. 
The  strains  e,  y in  Eq.  (l)'are  defined  with  respect  to  these 
local  cartesian  coordinates. 

The  geometry  of  the  shell  element  is  defined  by  approxi- 
mating the  global  cartesian  coordinates  of  a point  in  the  shell 
as  follows2,3 


(2) 


where 


X. 

1 


N. 


X 


global  coordinates  of  node  i on 
the  midsurface 

coordinate  normal  to  the  midsurface 
varying  from  -1  to  +1 
shell  thickness  at  nodal  points 
jth  component  of  nodal  unit  vector 
normal  to  the  midsurface 
interpolating  function  in  terms  of 
curvilinear  coordinates  £,  t]  embedded 
the  midsurface. 


in 
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Three  components  of  the  displacements  vector  in  the  global 


cartesian  coordinates  are  defined  a£ 


= E N . <V. 

1 ) 3J 


+ s I N.  t 
2 i 


all  a21 


12  ^22 


al3  a23 


<?1 

^2  X 


v;here  lb,  \b,  lb  = nodal  displacement  at  node  i 

.a^jf  ~ 3 component  of  the  nodal  unit  vectors 

a^  and  a2  tangent  to  the  midsurface 
<J>2  ” nodal  rotation  angles  of  the  normal 


-V 

vector  a. 


or  in  the  vector  form 


U = Aq+CBq  (4) 

where  q is  the  nodal  displacement  vector. 

To  obtain  strains,  first  we  need  the  derivatives  of  dis- 
placements U,  V,  W with  respect  to  n,  and  £ , and  the  inverse 


of  the  Jacobian 


J - 3 (x,y , z)/9 (C,n , C) 


For  example 


L 
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The  Jacobian  J is  dependent  on  £.  But  for  thin  shells  the  de- 
pendence of  J on  £ is  very  small  and  can  be  neglected.  As  an 
example  consider  a circular  ring  of  radius  R-  Ihen 


3X 

3C 


3Ni 

W 

3N 


*=  l 


3£ 


3C1(xi 

t-  a^.-,  1 
l 31^. 

(R  a31 

1 

t.  a0,  1 
i 31. 

r R a3i 

(l  + 

t. 

^ 2R^ 

(7) 


If  t:./2R  <<  1,  can  be  assumed  to  be  independent  of  £. 

In  this  case  the  strains  E in  the  global  cartesian  coordi- 

•*» 

nate  systems  can  be  expressed  as 


E = Ej  + 5 Ej  . (8) 

The  strains  e and  y in  the  local  coordinate  system  are  obtained 
from  E by  a proper  transformation.  The  direction  cosines 
v/hich  are  necessary  in  forming  the  transform  matrix  are  inter- 
polated from  the  components  of  a^,  a a^  vectois  at  the  nodes. 
Then  inolane  strains  e are  expressed  as 

)' 

> - Kj  q + C B?  q (9) 

The  dependence  of  Y on  t small  and  can  be  neglected  for 


thin  shells.  That  is 


-?3  3 


(10) 


For  the  present  element,  the  displacements  are  assumed  to 
be  biquadratic.  The  independent  transverse  shear  strains  are 
assumed  as  follows"1'^ 


y = a,  + a0  £ + a,  n 

'xz  1 2 3 

Y = a.  + ac  % + n 

ryz  4 5 6 


Yxz  = ct1  + a2  K + <*3  n + a4  C n 
V “ “5  + “6  e + a7  11  + °8  5 11 


(11) 


(12) 


VJritten  in  matrix  form 


y = p a 


(13) 


It  has  been  found  that  in  static  analysis  the  linear  strains 
in  Eq.  (11)  give  better  solutions  than  bilinear  strains  in 
Eq.  (12).^  Therefore  all  calculations  reported  in  the  pre- 
sent v/ork  were  performed  with  the  linear  transverse  shear 
strains. 
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Introducing  Eqs.  (4) , (9) , (10)  and  (13)  into  Eq.  (1)  , 
we  obtain, 

- l s{ST<5i  + 5i2  + *2>9 

+ aT  G q - j aT  H a}  (14) 

+ ^q^M^  + M^2  + M2)  q-<5qTQ=0 

where 

5i  » / b±t  9l  ?1  I J | dS  dn 

*12  “ ! «1T  Sl2  ?2  + ?2T  Sl2  5l>  lJl  « dn 

k2  = / b2t  C2  B2  |j|  dC  dn 
Si  - 1 9e 

Sl2  = f Se  l 

C2  = I Ce  £2  d? 

Mi  = / mi  AT  A | J | d£  dr) 

~12  = 1 ra12(~T  ! + ?T  |J|  dC  dn  (15) 

M2  = / m2  BT  B | J | d^  dn 
= / m d 5 

m12  = I 

m 2 = / m£2  d£ 

| J | = determinant  of  J 

G = / (PT  c b |J|  dc  dn 
~ An  ~ 

H = £ P c P J J | dC  dn 
n 

C3  = / Cy  d£  <5qT  Q = 6 W 


r 
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The  integration  through  the  thickness  is  carried  out  analyti- 
cally while  2x2  point  integration  is  used  in  £ - n plane. 
Taking  Sir  = 0 with  respect  to  a gives 


H a = G q 


a = H-1  G q 


Introducing  Eq.  (16)  into  Eq.  (14)  and  taking  6 tt  _ = 0 with 

mK 

respect  to  q leads  to 


E(Mi  + M12  + M2)q  + E(KX  + K12  + K2  + GTH-1G)q 
n ~ ~ n 


M q + K q = Q 


where  M = E (M^  + M^2  + M.>)  is  the  assembled  mass  matrix  and 
n T -1 

K = E (K^  + K12  + K2  + G H G)  is  the  assembled  stiffness  matrix 
n ~ 
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Fig.  C-l:  Coordinate  Systems  and  Nodal  Degrees 
of  freedom  for  a Shell  Element 


( 


59 


References 


1.  Proceedings,  "Third  Conference  on  Fibrous  Composites  in 
Flight  Vehicle  Design,"  NASA  TM  X-3377,  1976. 

2.  Griffin,  K.  E. , "An  Aeroelastic  Tailoring  Study  of  a High 
Aspect  Ratio  Wing,"  included  in  proceedings,  "Third  Con- 
ference on  Fibrous  Composites  in  Flight  Vehicle  Design," 
NASA  TM  X-3377. 

3.  Minich,  M.  D. , and  Chamis,  C.  C. , "Analytical  Displacements 
and  Vibrations  of  Cantilevered  Unsymmetric  Fiber  Composite 
Laminates,"  NASA  TM  X- 716 99,  1975. 

4.  Reissner,  E. , and  Stein,  M. , "Torsion  and  Transverse  Bend- 
ing of  Cantilever  Plates,"  NACA  TN  2369,  1951. 

5.  Ashton,  J.  E. , and  Whitney,  J.  M. , Theory  of  Laminated 
Plates,  Technomic  Publishing  Co.,  Stamford,  Conn.,  1970. 

6.  Leissa,  A.  W. , "Vibration  of  Plates,"  NASA  SP-160,  1969. 

^ 7.  MacNeal,  R.  H. , "A  Simple  Quadrilateral  Shell  Element," 

Computers  and  Structures,  Vol.  8,  pp.  175-183,  1978. 

8.  Spilker,  R.  L. , et  al. , "Static  and  Dynamic  Analysis  of 
Multilayer  Composite  Plates  and  Shells  by  the  Hybrid-Stress 
Finite  Element  Model,"  Aeroelastic  and  Structures  Research 
Lab  Report  ASRL  TR1S1-2,  M.I.T.,  June,  1976. 

9.  Lee,  S.  W. , "Finite  Element  Methods  for  Reduction  of  Con- 
straints and  Creep  Analysis,"  Ph.D.  Thesis,  Department  of 
Aeronautics  and  Astronautics,  M.I.T.,  February,  1978. 

10.  Bathe,  K.  J. , "Convergence  of  Subspace  Iteration,"  Paper 
presented  at  U.S.  - Germany  Symposium:  "Formulations  and 
Computational  Proceedures  in  Finite  Element  Analysis," 
M.I.T.,  August,  1976. 


1 


TABLE  1.  COMPARISON  OF  THE  FREQUENCIES  CALCULATED 
USING  MLP3K  AND  RS40  WITH  RITZ  SOLUTIONS 


HERE  B- BENDING,  T- TORS ION,  C-CHORDWISE 


TABLE  2.  FREQUENCIES  CALCULATED  USING  RS40 
FOR  A 3”  X 3"  x .0416"  PLATE 


TABLE  4.  CALCULATED  AND  OBSERVED  FREQUENCIES 
OF  ALUMINUM  PLATES 


* * 


PLATE 

MODE 

CALCULATED 

OBSERVED 

%DIFF 

SIZE 

(LENGTH  x 
WIDTH) 

FREQUENCY 

fc(Hz) 

FREQUENCY 

f0(Hz) 

3"  x 3" 

IB 

152.6 

153.0 

-0.3 

IT 

370.1 

367.1 

0.8 

2B 

937.9 

922.5 

1.7 

1C 

1193. 

1213. 

-1.7 

2T 

1356. 

1336. 

1.5 

6"  x 3" 

IB 

37.66 

37.  6 

0.2 

IT 

161.5 

158. 

2.2 

2B 

234.3 

234.9 

-0.3 

2T 

* 

523.6 

518.8 

0.9 

3B 

660.2 

658.1 

0.3 

FREQUENCIES  CORRECTED  TO  .0416" 


THICKNESS 
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TABLE  6. 

CALCULATED  AND  03SERVED  FREQUENCIES 

OF  6"  x 3"  8-PLY  GRAPHITE/EPOXY  PLATES 

LAMINATE 

MODS 

CALCULATED 
FREQ. (Hz) 

OBSERVED 

FREQ. (Hz) 

%DIFF. 

[02/+30]p 

IB 

(68.98) 

58.3 

(18.3) 

IT 

139.1 

148. 

-6.0 

2B 

408.7 

362.7 

12.7 

2T 

526.0 

508. 

3.5 

1C 

588.8 

546. 

7.8 

[0/+45/90 ] s 

IB 

(59.63) 

48.6 

(22.7) 

IT 

176.3 

169. 

4.3 

2B 

345.9 

303.0 

14.2 

2T 

592.2 

554. 

6.9 

1C 

820.4 

739. 

11.0 

[+45/+45 ] s 

IB 

(39.05) 

31.3 

(24.8) 

2B 

192.8 

185.8 

3.7 

IT 

229.3 

214. 

7.1 

3B 

565.8 

533. 

6.2 

2T 

707.6 

653. 

8.4 

NOTE:  VALUES  IN  PARENTHESIS  ARE  THOUGHT  TO  BE 

DUE  TO  COMPUTATIONAL  ERROR. 
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TABLE  7.  CALCULATED  AND  OBSERVED  FREQUENCIES 

OF  6"  X 3"  GRAPHITE/EPOXY/ALUMINUM  PLATES 


LAMINATE 

MODE 

CALCULATED 
FREQ. (Hz) 

OBSERVED 

FREQ. (Hz) 

%DIFF 

tVA^S 

IB 

113.8 

101.7 

11.9 

IT 

229.3 

229. 

0.1 

2B 

680.9 

631.9 

7.8 

2T 

885.8 

865. 

2.4 

1C 

1168. 

1129. 

3.5 

[ 0/+45/90/A1] * 

IB 

86.41 

75.9 

13.8 

IT 

314.6 

302. 

4.2 

2B 

506.1 

469.6 

7.8 

2T 

1035. 

983. 

5.3 

3B 

1418. 

1306. 

8.6 

[+45/+45/A1]* 

IB 

67.32 

58.3 

15.5 

2B 

355.9 

351.6 

1.2 

IT 

380.2 

358. 

6.2 

3B 

1029. 

1006. 

2.3 

2T 

1187. 

1113. 

6.6 

* THE  LAMINATE 

CONSIST  OF 

4-PLIES/. 040" 

ALUMINUM/ 4-PLIES 

1 


TABLE  8.  CALCULATED  AND  OBSERVED  FREQUENCIES  OF 


6"  x 3 
SHELL 

" GRAPHITE/EPOXY 
SECTIONS 

CYLINDRICAL 

LAMINATE 

MODE 

CALCULATED 
FREQ. (Hz) 

OBSERVED 
FREQ. (Hz) 

%DIFF. 

[02/+30]s 

IT 

169.9 

161. 

5.5 

IB 

291.6 

254.1 

14.8 

1C 

597.2 

555.6 

7.5 

2T 

719.5 

670.0 

7.  4 

2C? 

834.0 

794. 

5.0 

[0/+45/90 ] g 

IT 

182.2 

177. 

2.9 

IB 

231.2 

201.  8 

14.6 

2T 

703.6 

645. 

9.1 

1C 

809.5 

754. 

7.4 

? 

980.1 

884.8 

10.8 

[+45/+45]g 

IB 

144.7 

145.3 

-0.4 

IT 

248.3 

222. 

11.  8 

2T 

770.8 

712. 

8.3 

2B 

814.1 

774.2 

5.2 

1C 

1042.5 

997. 

4.6 

A 


( 


70. 


C 


FIG.  2 FLOW  DIAGRAM  FOR  PROGRAM  SHELL  VIBES 
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TG.  9 2nd  TORSION  (5th)  MODE  OF  AN 


ALUMINUM  PLATE 


CALCULATED  MODE  SHAPE  OBSERVED  MODE  SHAPE 

263.3  Hz  234.2  Hz 

FIG.  10  1st  BENDING  (1st)  MODE  OF  AN 
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CALCULATED  MODE  SHAPE 


225.4  Hz 


OBSERVED  MODE  SHAPE 
196.4  Hz 


FIG.  15  1st  BENDING  (1st)  MODE  OF  A 
[0/+45/90]  PLATE 
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140.8  Hz 

131. 

2 Hz 

FIG.  20  1st  BENDING  (1st) 

MODE  OF  A 

[+45/+451  PLATE 
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CALCULATED  MODE  SHAPE 
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OBSERVED  MODE  SHAPE 
790.5  Hz 


FIG.  22  2nd  BENDING  (3rd)  MODE  OF  A 
[+45/+45J  PLATE 
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FIG.  23  1st  CHORDWISE  (4th)  MODE  OF  A 
[+4S/+45]  PLATE 
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CALCULATED  MODE  SHAPE 
161.5  Hz 


OBSERVED  MODE  SHAPE 
158  Hz 


FIG.  26  1st  TORSION  (2nd)  MODE  OF  AN 
ALUMINUM  PLATE 
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FIG.  32  2nd  BENDING  (3rd)  MODE  OF  A 
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FIG.  34  1st  CHORDWISE  (5th)  MODE  OF  A 
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FIG.  35  1st  BENDING  (1st)  MODE  OF  A 
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FIG.  37  2nd  BENDING  (3rd)  MODE  OF  A 
[0/+45/90]  PLATE 
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FIG.  40  1st  BENDING  (1st)  MODE  OF  A 
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FIG.  53  2nd  TORSION  (4th)  MODE  OF  A 
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FIG.  55  1st  BENDING  (1st)  MODE  OF  A 
[+45/+45/A11  PLATE 
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FIG.  66  1st  BENDING  (2nd)  MODE  OF  A 
[Q/+45/90]  SHELL  SECTION 


V,V.  \ 

414  3)3 

f. r,  o f.  r. f.  r»  & r r.  *,  r,  < . ' ' r.  * . f.  r> 

1V.SS  1 

444  3 13 

SS55S  \ 

4444  ))) 

SSSVi  ! 

[1444  )).l 

sss 

SS'iSS 

14444  33: 

ww 

SSSr.V» 

v«4)i  y. 

ssssv. 

S5S 

5S5‘iSS 

3 4444  3. 

SSSWS'.S5 

S5SSSS5! 

mu 

s 

*iS^SSS5S 

SSV>S*V»S 

\ 1444 

V 

SSS^SSSSO',  sw 

S1SSS15S5 

V;  44  4 4 

I44l\ 

. SV.WWW5 

r*W.V.S5*S 

\ 4444 

444444 

X sv»ss‘>ss\,sr.*; 

WWWS5 

\ 4 m 

• 414)X  SSSV.'.SSSS 

OV^SWjSS 

34441! 

41  mA  SWWi 

WSSWISS 

1 

44444X  SV,SS'.V,SS!,Sr.r.S5S5 

\l< 

3)3 

44l4i\  55S5 

5Sr>,.,»'»55S55S5S  V 

3333 

44144\  SS'j*5SS5S5SSSSSS* 

34333 

iimX 

5SSSVi5S'.5Sr. 

SS5SS 

33331  4444434 

3333  44  44  4 X 

SiSISSWi 

3331  «*««\ 

S$S5SWSS5SSr.5S«5 

22 

3333  ••444\ 

SWSS5WSSSSWS‘ 

222 

3333  «4444\ 

SMWSSSWS1WV 

2222 

3333  44444'X 

S5WW 

2222 

3)3  4144)\ 

4SSS055 

SSV.SSf.W 

222 

3)3)  l<U4«A 

SS51M0 

S«55S‘liS50l 

2222 

3)3)  444441 

5S5S5SS5S5SS 

222 

3333  4«444>A 

SWSWW 

»n«nM 

444441 

I 

44144  I 

nm*  | 

4 14  lH 
<111)44 
44)144 
• 41444 
444141 
444444 


3)3)13 

414441  1 

31)3) 

4)444911 

CALCULATED  MODE  SHAPE 
703.6  Hz 


OBSERVED  MODE  SHAPE 
645.  Hz 


FIG.  67  2nd  TORSION  (3rd)  MODE  OF  A 
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FIG.  69  5th  MODE  OF  A [0/+45/90] 
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FIG.  72  2nd  TORSION  (3rd)  MODE  OF  A 
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FIG.  73  2nd  BENDING  (4th)  MODE  OF  A 


[+45/+451  SHELL  SECTION 
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FIG.  74  1st  CHORDWISE  (5th)  MODE  OF  A 
I+45/+45 ] SHELL  SECTION 


